Quantum transport: The link between standard approaches in superlattices 
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Theories describing electrical transport in semiconductor superlattices can essentially be divided 
in three disjoint categories: i) transport in a miniband; ii) hopping between Wannier-Stark ladders; 
and iii) sequential tunneling. We present a quantum transport model, based on nonequilibrium 
Green functions, which, in the appropriate limits, reproduces the three conventional theories, and 
describes the transport in the previously unaccessible region of the parameter space. 



OO 



i 

-a 
c 

o 
o 



> 
OO 

m 



OS 



i 

C 

o 
o 



Ever since the pioneering work of Esaki and Tsu |Q , 
which drew attention to the rich physics and potential 
device applications of semiconductor superlattices, these 
man-made structures have remained a topic of intense 
research. Semiconductor superlattices have proven to be 
a fruitful platform for studying a wide range of transport 
phenomena, such as their intrinsic negative differential 
conductivity 0, the formation of electric field domains 
||, Bloch oscillations Q, as well as dynamical localiza- 
tion jj| and absolute negative conductance H under ex- 
ternal irradiation, just to mention a few. 

These phenomena depend crucially on the relations of 
the energy scales involved, namely the zero-field mini- 
band width (which is four times the interwell coupling 
Ti), the scattering rate T/h, and the potential drop per 
period (= eFd, where F is the applied static field and d 
is the superlattice period). Three distinct approaches 
have been used to describe transport in the parame- 
ter space spanned by (Ti,eFd, T): miniband conduction 
(MBC) §||, Wannier-Stark hopping (WSH) §], and se- 
quential tunneling (ST) While the ranges of va- 
lidity of the different approaches have been addressed 
qualitatively before [[ll|--[l4[ , no explicit calculations have 
been presented where the different ranges can be iden- 
tified and the transition between them can be stud- 
ied. In the present Letter we present such a calculation, 
based on nonequilibrium Green functions. The calcu- 
lated current-field relations are shown to reflect the re- 
sults from the simple approaches (MBC, WSH, and ST, 
which will be reviewed below) in their respective ranges 
of validity sketched in the "phase-diagram" presented in 
Fig. [I. While similar diagrams have been obtained in 
Refs. ||@ from more phenomenological arguments, we 
will derive the borderlines from our Green function anal- 
ysis here. 

Now, we introduce the model assumptions which will 
be used in each of the following approaches. We restrict 
ourselves to the lowest miniband of the superlattice. Our 
basis set are orthonormal wave- functions ^ n (z)e l ^/A 
where the ^-direction denotes the growth direction. The 
^ n (z) = ^(z — nd) are localized in well n (for exam- 
ple one may use the Wannier- functions). Here r and k 
denote two-dimensional vectors within the (x, y)-plane 
(with area A) which is assumed to be separable from 



the z-direction. For parabolic dispersion Ek = fi. 2 k 2 /2m 
(with the effective mass m of the conduction band) we 
thus have a constant density of states po = m/nfi 2 per 
area and period. The single-particle part of the Hamil- 
tonian within nearest neighbor coupling is then given by 



H 
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^Ek-neFd). (1) 



Furthermore we consider a phenomenological scattering 
process at (5-potentials with density Nd and matrix el- 
ement 5 n-m V/A, leading to a scattering rate 1/to = 
To/h = NdTrV 2 po/h between the k-states within a given 
well. Finally, we assume that the in-scattering term is de- 
termined by a Fermi-distribution np(E) = [1 + exp((E — 
/i)//csT e )] _1 with electron temperature T e and chemi- 
cal potential p.. This assumption establishes internal en- 
ergy relaxation without specifying the corresponding pro- 
cesses in detail. It has been has been implicitly used in 
the standard approaches 10 as well. 



Miniband conduction (MBC): For zero electric field 
Eq. ([!]) is diagonalized by a set of Bloch functions <f q (z) = 
^2 n e mqd "$> n (z) and the dispersion relation is given by the 
miniband E(q) = 2TiCos(qd). The stationary Boltzmann 
equation for the distribution function /(g,k) is then 



eF df(q, k) _ n F (E(q) + E h ) - f(q, k) 



Ti dq 



T(E(q) + E k ) 



(2) 



where the relaxation-time approximation corresponds to 
our assumption on scattering mentioned above. For our 
scattering model, we obtain the relaxation time t(E) = 
t for E > 2\Ti\ and t(E) = ttt / arccos(-£y2|Ti|) for 
-2|Ti| < E < 2|Ti|. Eq. (g) is solved numerically and 
the current is calculated from 



J(F) 



d 2 k I dqf(q,k)- 
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^TT 3 h J J^^/ri * J ' dq 
The electron density per period is given by 



N 2D = — I d'k 



-TT/d 



dqf(q,k) 



(3) 
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and is used to determine the chemical potential (which is 
field dependent due to the energy dependence of t(E)) for 
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a given electron density. This approach can be extended 
beyond the relaxation time approximation |p^| , |l6| , but 
the generic features remain unchanged. 

Wannier- Stark hopping (WSH): In the presence of an 
electric field, the eigenstates of the Hamiltonian become 
the localized Wannier-Stark states, 



(5) 



with energy E v = —veFd where J n (z) is the Bessel func- 
tion of the first kind. Scattering causes hopping between 
the different states. Within Fermi's golden rule, the cur- 
rent is given by 



w = E<iE 
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d 2 fc [n F (E k ) - n F {E k + leFd)] . (6) 



Here the term E^[Jn 2 arises due to the spatial 
overlap of the Wannier-Stark functions and the Fermi- 
functions reflect our assumption regarding in-scattering. 
The electron density per period is given by: 



N 2 d = Pok B T e log 



I + exp 
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which relates fj, to N%d- Again, it is possible to general- 
ize this approach to more realistic scattering mechanisms 

Sequential tunneling (ST): In this approximation the 
phase information is lost after each tunneling event be- 
tween adjacent wells. The scattering within a well is 
treated self-consistently by solving for the spectral func- 
tions A(£,k); in this work we use the self-consistent 
Born-approximation pf | for the self-energy. The transi- 
tions to neighboring wells are calculated in lowest order 
of the coupling yielding 14 I9|j20[|: 



J{F) = ^J d2k J ^ T ? A (^k)A(£ + eFd,k) 

x [n F {£) - n F {£ + eFd)] . (8) 
The carrier density is given by: 

N2D= 2^J d2fc / £ ni?(£ )A{£ ' k) • (9) 

This approach gives quantitative agreement with exper- 
iments in weakly coupled structures when realistic mod- 
els for impurity and interface scattering are employed 

The important issue to recognize is that these three 
approaches treat scattering, external field, and coupling 
within different approximations. MBC does not properly 
include field-induced localization because of its inherent 
assumption of extended states, WSH treats scattering in 



lowest order perturbation theory (in particular, there is 
no broadening of the states), and ST is explicitly low- 
est order in the interwell coupling. In contrast to these 
shortcomings, a full quantum transport theory, based on 
nonequilibrium Green functions pl[ , is able to treat scat- 
tering, electric field, and coupling on equal footing. Such 
an approach was performed in Ref. p3[ using a basis 
of Wannier-Stark states and restricting the analysis to a 
high electron temperature. Here we work within the basis 
ty n {z) and consider the general situation which allows an 
analysis of transitions between the simplified approaches 
MBC, WSH, and ST. 

Nonequilibrium Green functions (NGF): Here the cur- 
rent and electron density are given by Jl9|j2l| ] 

J(F) = / d2 ^{ T i G »+i,»( f . t ,k)} (10) 

N 2D = ^ J d 2 kG< n (t,t,k) (11) 

where G< n (t, t', k) = i(oj 1 (f 1 k)o m (<,k)), and a£(i, k) 
and a n (t, k) are the creation and annihilation operators 
for the state \E , n (z)e^ k ' r ' ) /A in well n. We also need 
the retarded Green function: G^ n (t,t',k) = —iO(t — 
t')({a m (t, k), a}j(£', k)}), where {A, B} denotes the anti- 
commutator. In the stationary state the Green functions 
only depend on the time difference r = t — t' , and we 
define the Fourier transformation via pl| : 

G m , n (£,k)= J dre< £ - eFd ^)^G m , n (t,t-T,k). 

(12) 

Without scattering between the k-states and at T\ = 
the Green-functions are diagonal in the well index: 
Gmni^-i k) = i5m,ii3j, ct (^! k) with the free particle Green- 
function c£ et (£,k) = l/(£ - E k +i0 + ). The full Green 
function is then determined by the Dyson equation: 



G™(£, k) = g™ £ + eFd^— , k 



[S mtn + Y,KnA£ + eFd l —^,k 



KG \%[£ + eFd l -^,k) . (181 



Within the self-consistent Born approximation for the 
scattering the self-energy can be written as 



k) = <5 m ,„£; et (£, k) + T 1( Wm + Ti5, 



m—l,n 



(14) 

with £%*(£ ,k) = N d /A yj k , V 2 G^ n (£,k'). If the scat- 
tering term E^ ot is neglected, the solution corresponds to 
the Wannier-Stark states (||). On the other hand, ne- 
glecting the coupling T% gives the spectral functions used 
in the sequential tunneling model. Eqs. (0,n4h are solved 
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self-consistently for G Ict . Then G < is calculated via the 
Keldysh equation 

G< >n (S, k) = £ G£ mi ^ + eFd^p, k) 



X £< I £ + eFci mi 



■ f c : . r.,l ... m + n I k 



<G- V ,„(f + eFd^i^,kj ,15) 



According to our general assumption about in- 
scattering we replace (£ > k) by its equilibrium value 
— 2ini?(£)Im{E^ t (£, k)}. Finally, the current and elec- 
tron density are calculated via Eqs. (0,|ll|). The exten- 
sion of this model to more realistic scattering processes 
is straightforward by using the respective self-energies 
in Eqs. ( |l4| , |T5| ) and relaxing the assumption about in- 
scattering, although the calculations become very tedious 
(see, e.g., Ref. g2] where NGF has been applied to the 
resonant tunneling diode). 

Results: In Fig. || we display the evolution of the 
current-field relations for the different models from 
weakly to strongly-coupled superlattices. The curves for 
MBC, ST, and NGF are qualitatively similar for all cou- 
plings. For low electric fields the current increases lin- 
early with the electric field. Then there is a peak at 
intermediate fields and negative differential conductivity 
occurs at higher field. For small T\ as well as for high 
fields the result from ST is in quantitative agreement 
with the NGF result, while the results deviate for larger 
T\. In contrast, the result from MBC is in quantitative 
agreement with the NGF result for large T\ and small 
eFd. The WSH-result diverges for eFd -> @, but 
approaches the NGF result for large fields. These results 
as well as further calculations for ksT e < T,N 2D < Pa^ 
are summarized by Fig. [j], depicting the respective re- 
gions in parameter space, where the different approaches 
approximate the NGF result. 

Now we want to justify these ranges of validity by 
studying the quantum mechanical correlation between 
the wells n and m given by the retarded Green function. 
For a constant scattering self-energy S^ et (£,k) = —ir/2 
we have found an analytic solution of Eqs. (|lj,[l4|) 
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Here the terms G™±i n become of the order of G™^ when 
| J 1 « |J ±1 |, i.e., |4Ti/(eFd)sin(eFdr/2)| « . This 
can be used as an estimate for the boundary between 
localization and delocalization. Due to the exponential 
factor in Eq. ( |l7j ) only t < 2h/T is of relevance. If 
eFd > r the magnitude of the sine takes the average 
value w 1/V2. Then we find 2|Ti| » eFd. If, on the 
other hand, eFd < T we may replace sin(x) ~ x and 
have 2|Ti| w T at the time r = VZh/T. From these 
estimates we conclude that the states are essentially de- 
localized if 2|Ti| > T and 2|T\ | > eFd. In this case the 
miniband states form a useful basis as indicated in the 
upper left part of Fig. On the other hand for 2|Ti | < T 
or 2\Ti \ <C eFd the correlation functions G^* n vanish for 
to ^ n and the states are essentially localized so that the 
sequential limit can be used, as indicated in the lower 
part of Fig. |. 

For larger electron densities the agreement between 
the different approaches becomes better as shown in 
Fig. ||(a). These results together with further calcula- 
tions indicate that ST is also valid if N 2 d /pa ^ 2 ] Ti ] , and 
MBC is also valid if N 2D /po ^ T and N 2D /po > eFd. A 
similar trend is found for higher electron temperatures 
(Fig. H ( b)). This agrees with the analytic findings of 
Rcf. jn| where it is shown that NGF gives the same re- 
sult as MBC in the limit k B T e > |Ti|,eFd. 

In conclusion we have explicitly shown that a trans- 
port calculation based on nonequilibrium Green func- 
tions contains the simple approaches MBC, WSH, and 
ST as limiting cases. For low temperature and low elec- 
tron density the ranges of validity of the simplified ap- 
proaches are depicted in Fig. [I], while for higher electron 
densities or temperatures, these ranges are enlarged. The 
essential message of our analysis is that for wide regions 
in parameter space (but not everywhere!) a simplified 
theory can be found, which approximates the full theory 
satisfactorily. This should have important consequences 
for practical device modeling, where other complications, 
such as realistic scattering mechanisms, must be consid- 
ered as well. 

A.W. acknowledges financial support by the Deutsche 
Forschungsgemeinschaft . 



which is a superposition of broadened Wannier-Stark 
states. The Wannier-Stark ladder becomes resolved if 
eFd 3> T. This defines the region of validity for the 
WSH-approach, as indicated by the right region in Fig. ||. 
By Fourier transforming we obtain 



G^ n (M-T,k) = -ie(r)i 

X Jm — n 



i^ni+2L eF d-E k ) I e -rY/27i 



4Ti . (eFd 

sin r 

eFd V 2h 



(17) 
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FIG. 1. Regimes where the different transport models are 
valid for low electron densities and low temperatures. (For 
illustrative purpose we have translated the condition a ^> b 
in the text to a > 2b, where a and 6 denote the respective 
energy scales involved.) 




FIG. 2. Current-field relations calculated from nonequilib- 
rium Green functions (NGF) in comparison with the stan- 
dard approaches for N2D = 0.2Fopo and feT e = 0.2Fo. For 
Ti = 1.5r the current-field relation is also shown over a 
wider field range in the inset. Here one can see explicitly, 
that the NGF result leaves the MBC curve for eFd > Ti and 
approaches the ST and WSH curves for large fields. 
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FIG. 3. Current-field relations for (a) high electron density 
{N2D = 2ropo and ksTe = 0.2Fo) and (b) high electron tem- 
perature (fcsT e = 3Fo and N2D = 0.2ropo) for the different 
approaches. Note, that the MBC result deviates from the 
NGF result at eFd > N 2D /po in (a). 
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